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Imaging and spectroscopy with the Very Large Telescope have revealed 20 
Lya emitters within a projected distance of 1.3 Mpc and 600 km s~ x of the lu- 
1/-) ■ minous radio galaxy TN J 1338-1942 at z — 4.1. Compared to the field density 

of Lya emitters, this implies an overdensity on the order of 15. The structure 
has a projected size of at least 2.7 Mpc x 1.8 Mpc and a velocity dispersion of 
325 km s -1 , which makes it the most distant structure known. Using the galaxy 
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overdensity and assuming a bias parameter 6 = 3-5, the mass is estimated 
to be ~ 10 15 M Q . The radio galaxy itself is surrounded by an uniquely asym- 
metric Lya halo. Taken together with our previous data on PKS 1138-262 at 
z ~ 2.16, these results suggest that luminous radio sources are excellent trac- 
r£3 ■ ers of high density regions in the early Universe, which evolve into present-day 

clusters. The statistics of bright radio sources and of concentrations in the Ly- 
man break galaxy population are consistent with the picture that each of those 
concentrations harbours an active or passive luminous radio source. 
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Subject headings: galaxies: active - - galaxies: clusters: general - - galaxies: 
evolution — cosmology: observations — cosmology: early Universe 

1. Introduction 

Studies of high redshift (proto) clusters of galaxies (z > 2) can directly constrain theories 
of galaxy evolution and cosmological models (e.g. Bahcall & Fan 1998), but the detection of 
(proto) clusters at these redshifts using conventional optical and X-ray techniques is difficult. 
High redshift radio galaxies (HzRGs, z > 2) can help: they are among the most evolved 
and most massive galaxies in the early Universe (De Breuck et al. 2002) and are most likely 
located in dense (proto) cluster environments (e.g., Le Fevre et al. 1996; Pascarelle et al. 1996; 
Rottgering et al. 1996; Carilli et al. 1997; van Ojik et al. 1997; Pentericci et al. 2000b). In 
addition, HzRGs have properties that would be expected of forming central cluster galaxies. 
Their extremely clumpy morphologies revealed by Hubble Space Telescope images (Pentericci 
et al. 1999) are strikingly similar to simulations of forming brightest cluster galaxies, based 
on hierarchical models (e.g., Aragon-Salamanca, Baugh, & Kauffmann 1998). 

A pilot project on the Very Large Telescope (VLT) to search for Lya-emitting galaxies 
around the clumpy radio galaxy PKS 1138-262, resulted in the discovery of 14 galaxies and 
a QSO at approximately the same redshift as the radio galaxy (Kurk et al. 2000; Pentericci 
et al. 2000a). If the structure found is virialized, the total mass of the protocluster would 
be 10 14 M Q . Motivated by this result, we started a large program at the VLT to search for 
forming clusters (protoclusters) around HzRGs at redshifts 2 and higher. The first radio 
galaxy field we observed was TN J1338-1942 at a redshift of 4.1 (De Breuck et al. 1999). 
This HzRG is one of the brightest and most luminous in Lya known. Both its Lja profile 
and radio structure are very asymmetric (De Breuck et al. 1999), which indicates strong 
interaction with dense gas, and the rest-frame radio luminosity is comparable to that of the 
most luminous 3CR sources (Pusmhz — 4 x 10 35 erg s -1 Hz^ 1 sr _1 ). Here we report on the 
discovery of a substantial overdensity of Lja emitters around this radio galaxy at z ~ 4.1. 9 



9 Throughout this Letter, magnitudes are in the AB system and a A-dominated cosmology with H = 65 
km s _1 Mpc -1 , Qm = 0.3, and ft a = 0.7 is assumed. 
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2. Observations and candidate selection 

2.1. VLT imaging and selection of candidate Lya emitters 

We carried out narrow- and broadband imaging on 2001 March 25 and 26 with the 
8.2 m VLT Kueyen (UT2), using the imaging mode of the FOcal Reducer/low dispersion 
Spectrograph 2 (FORS2). At z — 4.10 the Lya emission line is redshifted to 6202 A, which 
falls in our custom narrowband filter with a central wavelength of 6195 A and FWHM of 60 
A. The broadband R filter had a central wavelength of 6550 A and FWHM of 1650 A. The 
detector was a SiTE CCD with 2048x2048 pixels, with a scale of / .'2 per pixel and a field 
of view of 6f8 x 6f8. We took 18 separate 1800 s exposures and one 900 s exposure in the 
narrow band and 21 exposures of 300 s in R, shifted by ~ 15" with respect to each other to 
minimize flat-fielding problems and to handle cosmic rays. The nights were photometric and 
the average seeing in both narrowband and broadband images was 0'.'65. The 1 a limiting 
magnitude per square arcsecond was 28.6 for the narrowband and 29.2 for the broadband 
image. For the flux calibration the spectrophotometric standard star GD 108 (Oke 1990) 
was used. The final images have sizes of 6'4 x 6'.2 (39.7 arcmin 2 ). The total volume probed 
at z = 4.1 by our narrowband filter is 7315 Mpc 3 . 

NTT imaging data, taken on 1998 April 27 and 29 under nonphotometric conditions 
with SuSI2, were also used to provide B- and /-band magnitudes for candidate emitters 
where possible. The 1 a limiting magnitude per square arcsecond was about 27.8 in both 
bands. 

For the detection and photometry of objects in the images, we used the object detection 
and classification program SExtractor (Bertin & Arnouts 1996). Detected objects had at 
least 9 connected pixels with values larger than the rms sky noise on the narrowband image. 
A photometric analysis was then carried out on both the narrowband and the broadband 
image. In total 2407 objects were extracted. Based on the statistics of our detection of 
Lja emitters around PKS 1138-262, we selected objects with a rest-frame equivalent width 
(EW ) greater than 15 A (or m B B— ^nb > 0.84) and significance £ > 3 as candidate Lya- 
emitting galaxies, with S the ratio of continuum-subtracted counts in the narrowband to the 
combined noise in the broadband and narrowband (Bunker et al. 1995). Of the 2407 objects, 
34 objects satisfied these criteria, including the radio galaxy. Of these 34 objects, 31 were 
detected in the i?-band image, 9 were detected in the NTT / image and 5 in the NTT B 
image. The 5 objects with a detection in B have a color of B — R ~ 0-1, inconsistent with 
the colors expected for galaxies at z = 4 (Steidel et al. 1999). Three of these objects are 
located in the halo of the radio galaxy, which affects the narrowband photometry. The other 
two are likely to be foreground objects with another line falling in the narrowband filter, e.g. 
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[O III] A5007 or [O II] A3727. Excluding the radio galaxy, the resulting 28 objects were our 
Lya-emitting candidates for follow-up spectroscopy. 

2.2. VLT spectroscopy 

On 2001 May 20, 21, and 22 we carried out spectroscopy using FORS2 in the mask 
multiobject spectroscopy mode with standard resolution. The nights were photometric with 
an average seeing of 1". The spectra were obtained with the 600RI grism with a dispersion 
of 1.32 A pixel -1 and a wavelength range from 5300 A to 8000 A. This grism was preferred 
because of its high throughput (peak efficiency is 87%). Two different masks were used to 
observe 23 of the candidate Lya emitters and the radio galaxy, with slit sizes of 1", resulting 
in a resolution of 6 A, which corresponds to 290 km s -1 at z — 4.1. The total exposure time 
was 31,500 s for the first mask and 35,100 s for the second mask. Wavelength calibration 
was obtained from exposures of He, Ar, Ne, and HgCd arc-lamps. The accuracy of the 
wavelength calibration was better than 0.05 A. For the flux calibration, long slit exposures 
with a 5" slit were used of the spectrophotometric standard stars EG 274 and LTT 7987 
(Stone & Baldwin 1983; Baldwin & Stone 1984). 

3. Results 
3.1. Line emitters in the field 

3.1.1. Line identifications 

Of the 23 candidates observed, 20 show an emission line with a peak between A = 
6187 A and A = 6216 A. The signal-to- noise of these lines is at least 10. Ten randomly 
selected spectra from our sample are shown in Figure 1. Two nondetections were very faint 
(m NB = 25.9 and m NB = 26.1 respectively). The third nondetection had a very low surface 
brightness. The success rate of our selection criteria was therefore 87%. 

The first question is whether the detected lines are indeed due to Lya at the redshift 
of the radio galaxy. Intervening neutral hydrogen will absorb emission blueward of the Lya 
line (Steidel et al. 1999). This discontinuity of the continuum over the Lya emission line is 
observed in one of the spectra, but the continuum emission of the other emitters was too 
faint to be detected (R ~ 27). Identification of the lines with [O III] A5007 at z ~ 0.24 can 
be excluded, because of the lack of confirming lines [O III] A4959 and FL3 in the spectra. 
The position of the emitter within the slit mask determines the wavelength coverage of the 
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resulting spectrum. In 9 cases this coverage was suitable to exclude the identification with 
[O II] A3727 at z ~ 0.66 on similar grounds. The detected emitters are anyway unlikely to 
be foreground [O II] galaxies. First of all, if one of the lines would be [O II], then the rest- 
frame equivalent width would be at least 70 A. A survey of nearby field galaxies, conducted 
by Jansen et al. (2000), gives a mean EWo of the [O II] line of ~ 30 A for galaxies with 
Mb = —16, which roughly compares to R ~ 27 at z = 0.66. Only 2 galaxies out of 159 
galaxies with [O II] in emission have EW ([O II]) > 70 A. The number of [O II] emitters 
expected in our field, using another study (Hogg et al. 1998), is ~ 7. Again, only a few 
percent of the [O II] emitters observed by Hogg et al. have an EW > 70 A. Therefore, from 
our sample of 20 emitters, < 1 is expected to be an [O II] emitter. 

Additional evidence that the observed lines are predominantly Lja lines at z ~ 4.1 
associated with the radio galaxy is provided by their velocity distribution, which has a 
dispersion of 326 ± 73 km s' 1 and a FWHM of 768 ± 172 km s" 1 (Fig. 2), which is 4 times 
smaller than the FWHM of the narrowband filter. Further, the velocity distribution peaks 
within 200 km s -1 of the radio galaxy Lya peak, corrected for absorption (De Breuck et al. 
1999). 

For all these reasons we interpret the observed emission lines as Lja. To estimate the 
redshifts, flux densities, and widths (FWHM) of the lines, a Gaussian function was fitted 
to each of the one-dimensional spectra. Details will be provided in a future paper (B.P. 
Venemans et al., in preparation). 

3.1.2. Significance and properties of the overdensity 

The next question to be addressed is to what extent the statistics of our detections 
represent a significant overdensity of galaxies in the field. A "blank-field" study of Lja 
emitters at approximately the same redshift is the Large- Area Lyman Alpha survey (LALA 
survey, Rhoads et al. 2000). Preliminary results of this survey indicate a number density of 
4000 ± 460 deg~ 2 Az^ 1 for objects with EW > 15A and line + continuum > 2.6 x 10~ 17 
erg s _1 cm~ 2 . For our field the second criterion corresponds to m^B < 24.55 and we expect 
2.3 ±0.3 such Lya emitters within our probed volume. However, 9 of the confirmed emitters 
satisfy the above criteria (including the radio galaxy). In our cosmology the comoving volume 
density of the LALA survey is ulala = (3.1 ±0.4) x 10~ 4 Mpc -3 . The volume density in our 
field is nigss = 9/7315 Mpc" 3 = (12 ±4) x 10~ 4 Mpc~ 3 . The difference in number density is 
^1338 I^lala = 4.0 ± 1.4. However, the FWHM of the velocity distribution is approximately 
4 times smaller than the FWHM of the filter (see Fig. 2). This implies that our radio galaxy 
field is overdense in Lyct emitters by a factor of 15 compared with a blank field. 
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The spatial distribution of the emitters is not homogeneous (Fig. 3). The structure 
appears to be bound in the northwest but our FOV is not large enough to show a boundary 
in the south. The size of the structure is therefore at least 6'x 4', which corresponds to greater 
than 2.7 Mpc x 1.8 Mpc. Remarkably, the radio galaxy probably the most massive system 
in the structure, does not appear to be at the center of the protocluster, in disagreement 
with models for the formation of dominant cluster galaxies (e.g., West 1994). 

3.2. Radio galaxy halo 

The radio galaxy is located close to the apparent north-west boundary of the galaxy 
overdensity structure. The radio emission is dominated by two components, separated by 5'.'5 
(De Breuck et al. 1999). The brightest component coincides with the optical emission, while 
the other is in the south-east, towards the center of the galaxy overdensity. A spectacular 
feature of the radio galaxy, visible in the narrowband image, is the large Lja halo (Fig. 
4). Although giant Lya halos are a common feature of HzRGs, they are usually fairly 
symmetrically extended around the radio galaxy. In the case of TN J1338-1942, the halo is 
highly asymmetric and extends for ~ 15" (110 kpc) to the north-west, i.e., away from the 
center of the overdensity structure. In its narrowness and asymmetry the TN J1338-1942 
halo is unique among known distant radio galaxies. Possible mechanisms to produce this 
structure include cooling flows in colliding sub-structures and buoyancy effects (Gisler 1976) 
and will be considered in a subsequent paper. 

4. Discussion 

4.1. Nature of the overdensity 

Could the structure that we have detected be a protocluster at z ~ 4.1 that will evolve 
into a rich cluster of galaxies in the local Universe? At z — 4.1 the Universe is ~ 1.6 Gyr 
old, too short for the structure to have virialized since the mean crossing time for galaxies at 
the observed velocity dispersion is at least 4 Gyr. Thus the observed velocities are probably 
infall velocities of the galaxies accreting onto a large overdensity. The total mass of this 
structure can be estimated by using M = pV(l + 5 m ) with p the mean density of the 
Universe and S m the mass overdensity within our volume V (Steidel et al. 1998). The mass 
overdensity is related to the galaxy overdensity 5 g3i \ through 1 + b5 m — C(l + 5 ga i), where 
C takes into account the redshift space distortions caused by peculiar velocities and b is 
the bias parameter. From the statistics of redshift "spikes", Steidel et al. (1998) argue that 
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b > 4. Taking b in the range 3-5 and assuming that the structure is just breaking away from 
the Hubble expansion, 5 m is estimated to be 0.5-2.3, giving a mass of our structure of (1-2) 
xlO 15 M Q . This is comparable to the mass of the Lyman break galaxy overdensity (spike) 
found by Steidel et al. and to that of the Coma cluster (e.g., Fusco-Femiano & Hughes 1994). 

4.2. Relation to overdensity spikes 

It is instructive to compare the number of luminous radio sources at z ~ 3 to the number 
of redshift spikes. Steidel et al. (1998) estimate that about one velocity spike is detected in 
each 9'x 18' field in their spectroscopic survey, corresponding to 9 x 10 5 redshift spikes in 
the whole sky. Their survey is sensitive to redshifts between 2.7 and 3.4, corresponding to a 
cosmic evolution time of 0.6 Gyr. 

How many luminous radio sources are there at z ~ 3? Using the pure luminosity 
evolution model of Dunlop & Peacock (1990) to describe the steep spectrum radio luminosity 
function, we estimate that in the redshift range 2.7 < z < 3.4 there should be ~ 1.2 x 10 4 
radio sources with luminosities exceeding 10 33 erg s" 1 Hz -1 sr _1 at 2.7 GHz ("Cygnus A 
type" radio sources: P 2. 7gh z ( Cygnus A) ~ 2 x 10 33 erg s -1 Hz -1 sr _1 , Becker, White, & 
Edwards 1991). Assuming that HzRGs are only once active for 10 7 yr (Blundell & Rawlings 
1999), we expect the number of (previously) active radio sources in this redshift range to be 
~ 7 x 10 5 . 

Hence, the luminosity functions and lifetimes of luminous radio sources are consistent 
with every velocity "spike" in the space densities of Lyman break galaxies being associated 
with a massive galaxy that has been or will become a luminous radio source once. Note 
that West (1994) presented similar statistical evidence to argue that distant powerful radio 
galaxies are the precursors of cD galaxies at the centers of galaxy clusters. 

5. Conclusion 

We have found a structure of 20 Lya-emitting galaxies around the high redshift radio 
galaxy TN J1338-1942. The overdensity of this protocluster is on the order of 15 compared 
to field samples. Our results demonstrate that by z — 4.1, megaparsec-scale structure had 
already formed. 

Together with our previous data, this implies that the most luminous radio sources are 
tracers of regions of galaxy overdensity in the early Universe. The estimated masses of 10 14 
- 10 15 M & are consistent with the overdensities being ancestors of rich clusters of galaxies 
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in the local Universe. 
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Fig. 1. — Spectra of 10 of the 20 confirmed Lya emitters. For clarity, each spectrum is offset 
by multiples of 2.5 x 10~ 18 erg s _1 cm~ 2 A -1 . 
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Fig. 2. — Velocity distribution of the Lya emitters. The bin size is 175 km s _1 and the 
median of the velocities is taken as zero point. The velocity of the radio galaxy, corrected for 
absorption (De Breuck et al. 1999), is indicated by an arrow. The normalized transmission 
curve of the narrowband filter is also plotted. Note that the velocity distribution of the 
detected emitters is substantially narrower than the filter width and centered within 200 km 
s _1 of the redshift of the radio galaxy. 
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Fig. 3. — Spatial distribution of the 20 confirmed Lja emitters at z ~ 4.1 (circles) and the 
radio galaxy (square). The size of the circles is scaled according to the Lya flux of the object, 
ranging between 0.3 and 4.1 xl0~ 17 erg s -1 cm -2 . The structure appears to be bound in 
the northwest of the image and unbound in the south. Note that the radio galaxy is not 
centered in the galaxy distribution. 
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Fig. 4. — Continuum-subtracted Lya image of the radio galaxy halo. The contours represent 
the Lya flux density in the center of the halo. The surface brightness ranges from ~ 0.07 
to 1 x 10~ 17 erg s _1 cm" 2 arcsec -2 . The low-brightness halo is extended to the northwest, 
pointing away from the overdensity structure. This Lya halo is the most asymmetric radio 
galaxy halo known. 



